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We have extended the application of quantum renormalisation group (QRG) to the anisotropic 
Heisenberg model with next-nearest neighbour (n-n-n) interaction. The second order correction 
has to be taken into account to get a self similar renormalized Hamiltonian in the presence of n-n- 
n-interaction. We have obtained the phase diagram of this model which consists of three different 
phases, i.e, spin-fluid, dimerised and Neel types which merge at the tri-critical point. The anisotropy 
of the n-n-n-term changes the phase diagram significantly. It has a dominant role in the Neel-dimer 
phase boundary. The staggered magnetisation as an order parameter defines the border between 
fiuid-Neel and Neel-dimer phases. The improvement of the second order RG corrections on the 
ground state energy of the Heisenberg model is presented. Moreover, the application of second 
order QRG on the spin lattice model has been discussed generally. Our scheme shows that higher 
order corrections lead to an effective Hamiltonian with infinite range of interactions. 
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I. INTRODUCTION 

Quantum phase transition has been one of the most interesting topics in the area of strongly correlated systems in 
the last decade. It is a phase transition at zero temperature where the quantum fluctuations play the dominant role 
0. Suppression of the thermal fluctuations at zero temperature introduces the ground state as the representative of 
the system. The properties of the ground state may be changed drastically shown as a non-analytic behaviour of a 
physical quantity by reaching the quantum critical point. This can be done by tuning a parameter in the Hamiltonian, 
for instance the magnetic field or the amount of disorder. The study of the ground state and its energy is thus of 
central importance for understanding the critical behaviour of such systems. 

The technique of renormalisation group (RG) has been so devised to deal with these multi-scale problems 0, ■ 
In the momentum space RG which is suitable for studying the continuous systems, one itcrativcly integrates out small 
scale fluctuations and renormalizes the Hamiltonian. In the real space RG, which is usually performed on the lattice 
systems with discrete variables (i.e quantum spin chain) , an original Hamiltonian is replaced with an effective one for a 
lower energy subspace iteratively. In this approach the Hamiltonian is divided into inter-block (H BB ) and intra-block 
parts {H B ), H B is diagonalized exactly and then H BB is projected into the low energy subspace of H B 0. The 
accuracy of this method is determined by the number of states kept in the H B subspace and the approach to consider 
the effect of neglected subspace. The Ising model in a transverse field @ and the anisotropic Heisenberg model 7] have 
been studied by quantum renormalisation group (QRG) approach which gives the correct phase diagram. Moreover, 
the recent study on a more general model, XYZ in a transverse field, supports the power of this method to study the 
collective behaviour of the spin models yj. 

In this paper we are going to study the effect of higher order corrections on the QRG scheme. In this respect we will 
consider the one dimensional S = h antiferromagnetic XX Z chain with next-nearest neighbour (n-n-n) interactions. 
Because, even in the case of nearest neighbour Heisenberg model the n-n-n interaction will be generated for the 
renormalized Hamiltonian if we add the second order corrections. Moreover, we will obtain the phase diagram of this 
model which is a function of the anisotropics and the n-n-n interactions. We have calculated the effective Hamiltonian 
up to second order in H BB . The second order correction improves the accuracy of the results. However, it must 
be taken into account to get a self similar Hamiltonian after each step of QRG for the n-n-n XX Z chain. In this 
approach we have considered the effect of the whole states of the block Hamiltonian which are partially ignored in 
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the first order approach. The present scheme allows us to have the analytic RG equations, which gives a better 
understanding of the behaviour of system by running of the coupling constants. We have succeeded to obtain the 
phase diagram which contains the critical surface between the spin-fluid, Neel and dimer phases. The boundaries 
between these phases merge at the tri-critical point. The projection of our phase diagram on the A = 5 plane (the 
same anisotropy for nearest and next-nearest neighbour interactions) is in good agreement with the numerical ones 
. The QRG equations show that the anisotropy in the n-n-n term changes the phase boundaries significantly. We 
have also derived the staggered magnetisation in the z direction, that is the proper order parameter to show the phase 
transition between Neel-dimer and fluid-Neel phases. 

In the recent years, several interesting quasi-one-dimensional magnetic systems have been studied experimentally 
[lOl ITTl IT^ I . Among them, some compounds containing CuO chains with edge-sharing CuO^ plaquette were expected 
to be described by the XX Z model with next-nearest neighbour interaction The Hamiltonian of this model on 
a periodic chain of N sites is: 

J N 

H = 4 { £[°f<u + °M+i + A + + °W + 2 + s « +2 )]}, (i) 

i=l 

where J > and J x J 2 > are the first and second-nearest neighbour exchange couplings and the corresponding 
easy-axis anisotropics are defined by A and 8. For J 2 — 0, the ground state properties are well known by the Bcthe 
ansatz and the bosonization technique. In this case, for A < —1, the system is in the ferromagnetic phase, while 
for A > 1, it enters the Neel phase where the twofold degenerate ground states are separated from the excited ones 
by a finite gap ^{|. In the case of AY"-type anisotropy (—1 < A < 1) quantum fluctuations destroy the long range 
order even at the zero temperature. The ground state is characterised by gapless excitations and algebraic decay 
of spin correlations (spin-fluid) . A kind of frustration can be introduced to this model by adding the next-nearest 
neighbour interaction (J2 7^ 0) as well as the nearest neighbour ones. On the J2 = ^,A = <5=1 line, the model is 
known as the Majumdar-Ghosh Hamiltonian where the exact ground state has been obtained to be purely dimerised 
[l6j| . The dimer state is characterised by the excitation gap, the exponential decay of the spin correlation functions 
and the dimer long-range order. According to the above facts, there should exist the fluid, dimer and Neel phases 
in this model. The isotropic model have been studied intensively by analytic and numerical approaches |l6 | -|26 | . 
However, the QRG approach to this model has not been considered yet which gives a clear phase diagram in the 
presence of anisotropics. Moreover, as mentioned earlier, we are going to address the implementation of higher order 
RG corrections and discuss its features. We will explain the QRG scheme in the next section where the second order 
effective Hamiltonian and the renormalisation of the coupling constants are obtained. In Sec. Ill , We will present 
the phase diagram and explain the border between different phases. Finally, we will discuss on the features of the 
second order approach and higher order corrections. 



II. RG EQUATIONS 

The main idea of the RG method is the mode elimination or thinning of the degrees of freedom followed by an 
iteration which reduces the number of variables step by step until reaching a fixed point. We have implemented the 
Kadanoff 's block method for this purpose, because they are well suited to perform analytical calculations in the lattice 
models and they are conceptually easy to be extended to the higher dimensions. In Kadanoff's method, the lattice is 
divided into N' blocks with n s sites each and the original non-diagonal Hamiltonian is replaced with an effective one 
for a lower energy subspace. The Hamiltonian can be written 

H = H B + XH BB , (2) 

where the block Hamiltonian H B is a sum of commuting Hamiltonians, each acting on every block and A is a coupling 
constant which is already in H or else it can be introduced as a parameter characterising the inter-block coupling 
which can be set to 1 at the end of calculation. An exact implementation of this method j5j is given by the following 
equation, 

H e "=P B HP B , (3) 

where Pb is the projection operator. Eq.© suggests that we should search for the solution of the Eq.J2J| in the form 
of a perturbative expansion in the inter-block coupling parameter A, namely 

Pb = Pa + APi + X 2 P 2 + • • • , H eff = H^ f + \H{ ff + \ 2 H e 2 if + ■■■ . (4) 
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FIG. 1: The decomposition of chain into the three site blocks Hamiltonian (H B ) and the inter-block Hamiltonian (H BB ) 



To the zeroth order in A Eq.J2J) becomes 



H e Q n = P H B P . 



(5) 



Since H B is a sum of disconnected block Hamiltonians 



N' 



H B =Y,hf 



1=1 



one can search for a solution of Pq in a factorised form 



TV' 
1=1 



It can be found that 



»=1 

where (i = 1, • • • , k) are the k lowest energy states of hf. First and second order corrections are obtained to be 

Hf f = P Q H BB P , H e 2 ff = P Q H BB P U (6) 
where Pi is defined by the following equation 

Pl = {l- Po) „ * (1 - Po)H BB P Q . 



E -H B 

The second order correction can be written p7l |2S|. 



ff 2 e// = P [ J ff BS (l - P ») Eo 1 h b ( 1 - Po)H BB ]Po- 



(7) 
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We have considered three-site block (Fig.JIJ) and kept the degenerate ground states (\i>o), \ip'o)) of each block to 
construct the projection operator (P = \ipo)(ipo\ + IV'oKV'oD- The ground states and their corresponding eigenvalues 
of each block are 



|V*> 
Wo) 
eo ; 



=(| tU> + 9| TIT) 



1 



J 



=d ru> + 9i in) 



ITT)), 
IIT)), 



[J 2 - A - V(^2 



A 2 ) 2 



where A 2 = J 2 <5. In this notation | f), II) are the a z eigenbases, and 



q = --(J a + A - A 2 + V(^2 + A- A 2 ) 2 + 8). 

The interaction between blocks defines the effective interaction of the renormalised chain where each block is considered 
as a new single site. Calculating the effective Hamiltonian to the first order correction leads to the XX Z chain without 
next nearest neighbour interaction ( J' 2 = 0), i.e the effective Hamiltonian is not exactly similar to the initial one. The 
next-nearest neighbour interaction is the result of the second order correction. When the second order correction is 
added to the effective Hamiltonian, the renormalised Hamiltonian apart from an additive constant (cb) is similar to 
Eq.Q with the renormalised couplings. Thus, the effective Hamiltonian including the second order correction is 



H eff = _ 



J' 



N/3 



N/3 



N 



The renormalised coupling constants are functions of the original ones which are given by the following equations 

2 o > J 2 , A 



J' 



A: 



-Y( g ' + 2J i q)- T { 



)( 



f 



I 2 

— ( 

4 e - ei 



)( 



1 



e - e 2 ' ' 2 + q 2 

2, 



(8) 



) (p+q) x 



J 



(2 + q 2 )(2+p 2 ) 
( pq )( P + q + 4J 2 )( ( A - 2A 2 )pq + 4A 2 ) 
4 w 1 g 



)( 



) a ( 



4 2e - ei - e 2 2 + ^ 2 + p 2 ' 
( P + q + 2J 2 )((A-A 2 )pq + 2A 2 



A' = J( 



^) 2 (Ag 2 + 2A 2 (2- g 2 ) , 



(9) 



I 2 1 
— ( 



(p + g ) 2 + 4J 2 (p + g ) 2 

A f o i _2 \C o i „2 \ / 



4 v eo-ei /v ( 2 + g 2 )( 2+p 2 ) 
9 2 s2 j2 / 1 n/ l + 2J 2 g 



4 ( eo -e 2 )( 2(2 + g 2 ) )2 4 ( e - e3 )( 2 + g 2 )2 



J 



( 



4 2eo - ei — e 2 2 + p : 
2 



)( _I_ )( g ( P + g + 2J 2 ) )2 



2 + g 2 
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qj i + g ^2 ) , 2 

4 v 2e - e 2 - e 3 /v 2 + g 2 



)( 
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, J 2 {pq + q 2 + 2)+p + q 2 . 



4 v 2e -ei-e 3 /v (2 + g 2 )(2 + p 2 )V2 
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4 = 



j 2 

— ( - 

4 v 2 



( J 2 ( Sq + p ) + q(p + q ) ) 5 



( J 2 ( 1 + g 2 ) + g) 2 

eo - e 3 



(eo-ei )( 2 

( g 2 + J 2 g ) 2 
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+ P 2 ) 



(10) 



/J'. 



A 2 = 



/2 

— ( - 

4 v 2 



( A 2 q(p-pq 2 +q ) + f W 3 ) 2 
(co-ei )( 2+p 2 ) 



(11) 



( Ag 2 + A 2 ( 2-g 2 
2( e - e 2 ) 



/J'. 



in which 



en = je (12) 
+i( „( A + 2A S ( 2 - „ ) ) f ) + ( -L- )( 5(T L FT f { f + 8 A* ) 



_J_,<^t«!,> + ( » xjh^jh 

e - e3 2 + 2eo - ei - e 2 2 + ^ 2 + 



( ( (p + q)q + 2J 2 qY + 2(Apq + A 2 {2-pq) f ) 

"(o ; )(^) 2 (^)((p + ?) + ^(p + «z) + 2J 2 ) 2 

2e - ei — e 3 2 + q z 2 + p A 



+ ( 5 " )(^) 2 ( 1 + ^ ) 2 

2e - e 2 - e 3 2 + q- 



3 16 ^ eo - e 3 2 + 
1 



v )( ~ ) 3 ( 2 ( « 2 + Jil f + ( A q 2 + 2A 2 - A 2 q 2 f 

e - e 2 2 + <^ 

( )( ) 3 ( T ^)(&[q 2 +P Q + q + J2 P} 2 

eo — ei 2 + V 

[p?(A g 2 + 2A 2 - A 2 g 2 ) + A 2 q 2 (2-pq) ] 2 ) }. 



P = J 2 + A - A 2 - yJ(J 2 + A-A 2 ) 2 



In the above equations, ei, e 2 and e 3 are the first, second and third eigenenergies of the block's Hamiltonian which 
have analytic expressions in terms of the coupling constants but not presented here. Due to the level crossing which 
occurs for the eigenstates of the block Hamiltonian, RG equations are valid for i [(A — A 2 ) — y(A — A 2 ) 2 + 4] < J 2 < 
i[(A — A 2 ) + y/{A — A 2 ) 2 + 4]. We have plotted the RG flow and different phases in Fig.Q which will be discussed 
in the next section. 



III. PHASE DIAGRAM 



The RG equations show the running of J coupling to zero which represents the renormalisation of energy scale. At 
the zero temperature, phase transition occurs upon variation of the parameters in the Hamiltonian. In the region of 
planar anisotropy < A < 1, the nearest neighbour interaction (J 2 = 0) is known not to support any kind of long 
range order and the ground state is the so called spin-fluid state. Increasing the amount of anisotropy is necessary to 
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FIG. 2: Phase diagram of the XXZ model with next-nearest neighbour interaction and different anisotropies. Open circles 
show the fixed points, the tri-critical point (A| = J| — 0.155, A* = 1) and the XY fixed point (A2 = J2 = A* = 0) . For 
< A < 1, the front side of diagram is the spin-fluid phase (I) and the back side is the dimer phase (III). The checkerboard 
pattern represents the boundary between them. In the case of A > 1 there are two phases, the Neel-phase (II) which exists 
below the semi-circular lines and the dimer-phase (III) which is above them. The boundary between Neel and dimer phases is 
shown by the semi-circular lines. The surface which is depicted by the dash-dotted lines is the boundary between the dimer-Neel 
phases which is a vertical plane at A = 1. The phase transition between the spin-fluid phase (I) and the dimer phase (III) 
occurs at the vertical plane A = 1 which has not been plotted in this diagram to avoid it being complex. 



stabilize the spin alignment. For A > 1 the ground state is the Neel ordered state. In the case of J2 > 0, the nearest 
neighbour and next-nearest neighbour couplings are in competition with each other. The latter thus frustrates the 
ordering tendency of the former. For < A < 1 the interplay of the two competing terms in the presence of quantum 
fluctuations produces the dimer phase for J2 > Jf (A, A2). Our RG equations show that the phase boundary between 
the dimer and the spin-fluid phase depends on the n-n-n anisotropy coupling (A 2 ) and can be described by a two 
dimensional surface which is convex in a view from the spin-fluid phase (see Fig. © , the checkerboard curved plane) . 

The dimer or spin Peierls phase has a spin gap and a broken translation symmetry (the unit cell is doubled) in 
the thermodynamic limit. The dimcr-fluid transition is known to be of Bcrzinskii-Kosterlitz-Thouless (BKT) type 
|2flj . At the BKT transition point, the divergence of the correlation length is not of the usual power law type but 
very singular 30]. In fact, all the derivatives of the inverse of the correlation length at the critical point are zero. 
Moreover, at the BKT critical point there appear the logarithmic correction as a finite size effect which converges 
very slowly, in various quantities, such as correlation function and susceptibilities, therefore it is very difficult to find 
the critical point of the BKT type transition accurately. Two effective methods have been introduced to find this 
critical point. In the work of Roomany et. al |3lJ, (3- function is calculated numerically, while in the Nomura-Okamoto 
approach the fluid-dimer transition is determined by the degeneracy between the doublet excitation and the dimer 
excitation . Moreover, in Ref.[26| a proper structure factor has been introduced to probe this transition. However, we 
determine the fluid-dimer phase transition by using the running of couplings under RG. In the spin- fluid phase, the 
anisotropy and next-nearest neighbour couplings are irrelevant while in the dimer phase they run to the tri-critical 
point (Af = J| — 0.155, A* = 1). For small A2 the fluid-dimer phase boundary (Jf (A)) shows an inclination to the 
lower values versus A while it behaves conversely for higher A2 as shown in Fig. © • The comparison with numerical 
results 9] for A = 6 shows good qualitative agreement. The projection of the three dimensional phase diagram 
(Fig.©) on the A2 = J2A plane is shown in Fig.©. Different phases and the running of RG flows can be observed 
simply. However, from quantitative point of view we got higher values for Jf than the numerical results. For instance, 
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FIG. 3: The projection of the 3-dimensional phase diagram (Fig.JJ) on the A2 = J2A plane. Arrows show the running of 
couplings under RG. Open circles show fixed points. The dashed lines in the spin-fluid and Neel phases are the lines of the 
Gaussian fixed point. Different phases are labeled by: (I) spin-fluid, (II) Neel order, (III) dimer. The Filled circle denoted by 
J2 = 1/2 on the A = 1 line is the pure dimer state. 



at A = the RG analysis gives J| — 0.44 which can be compared with the numerical result of ~ 0.33 presented in 
Ref.0- The difference is inherited to the QRG scheme where the boundary condition of the isolated block does not 



In the region of the antiferromagnetic anisotropy (A > 1) where the Neel phase exists, it is destabilized by sufficiently 
strong competing next-nearest neighbour couplings, J| and AJj. The Neel phase appears just by crossing the A = 1 
plane at A2 = and J2 = 0. In the A2 = plane and for A > 1 the model will pass through a phase transition from 
Neel to dimer phase for J2 > J|(A). This boundary is the intersection of the semi-circular lines and the A2 = plane 
which happens around J2 ~ 1.5 (see Fig. 101 ). The Neel ordered is also broken by increasing the anisotropy of the 
n-n-n interaction. This is a phase transition to the dimer phase which is shown in Fig. (O by the semi-circular lines 
where the dimer phase exists above them. Thus, both the n-n-n exchange (J 2) and n-n-n anisotropy (A2) drive the 
model from the Neel ordered phase to the dimer one. In this respect, the boundary between the Neel and the dimer 
phase looks like an arcade in the phase diagram. 

We have probed the boundary of the Neel-dimer transition by calculating the staggered magnetisation in the 
^-direction as an order parameter (Fig.Q), 



The staggered magnetisation (Sm) is zero in the dimer phase and has a nonzero value in the Neel phase. Thus the 
staggered magnetisation is the proper order parameter to represent the Neel-dimer transition. We have plotted Sm 
versus A 2 for different values of A = 1.05, 1.15 and J2 = 0.3, 0.6 in Fig.QJ. The staggered magnetization goes to zero 
continuously at a critical value of A^A, J2) which shows the destruction of Neel order. 

A similar calculation has been done to trace the phase transition between the Neel and the spin-fluid phase. The 
crossing boundary of the fluid and Neel phases always stays at A = 1. The reason is related to the SU{2) symmetry 
of the Hamiltonian on the line A = 5 = 1. The QRG approach preserves this symmetry which can be seen by the 
irrelevant direction of A2 = J 2 in the A = 1 plane toward the tri-critical point. 

The significant result of our calculations occurs at the isotropic plane, A = 1. For small J2 the RG equations show 
running of J2 to zero except at the isotropic plane (A = 1). This means, if we start with the XXX model ( J2 = 0), 
the next-nearest neighbour interaction (J 2 ) is generated and runs to the tri-critical point, Aj = J| = 0.155, A* = 1. 
It might be interesting to mention that this fixed point corresponds to the reported experimental value for CuGeO^ 
(J2 = 0.13) It is worth mentioning that for the isotropic case A = 5 = 1, if the projection operator for the 
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FIG. 4: The staggered Magnetisation in z-direction for different values of A and J2 versus A2. The order parameter (Sm) goes 
continuously to zero which represents the transition from Neel to dimer phase. 



second order RG, (1 — Pq), does not contain all of the remaining states of the block we get J| = 0.0761 |2jj. This 
shows the increase of J r | when we have considered the correlation effects more precisely. Thus, we expect a nonzero 
value for J 2 if we go further to consider the higher order RG approximation. This is also true for J 2 = as the initial 
value, the corresponding fixed point is nonzero. So, the manageable fixed point, the fixed point which enables us to 
get the physical quantity, represents the next-nearest neighbour interaction for the isotropic Heisenberg chain. At 
this fixed point the correlation length is zero and the system behaves as a classical one. It represents a frustrated 
chain of classical spins which does not posses any ordering. It might be an explanation why the nearest neighbour 
one-dimensional antiferromagnetic Heisenberg chain does not show a true long range order. 

The other remarkable result of our RG flow is the observation of the specific lines in the spin-fluid and Neel phases 
which has been shown by dashed lines in Fig.®. We have linearised the RG flow at (A 2 = J 2 — 0.155, A* = 1), 
and found one relevant and two irrelevant directions. The eigenvalues of the matrix of linearised flow are Ai = 1.24, 
A 2 = 0.46 and A3 = 0.41. The corresponding eigenvectors in the |A, J 2 , A 2 ) coordinates are |Ai) = | — 0.98, 0.02, —0.20), 
|A 2 ) = I - 0.88,0.15, -0.44) and |A 3 ) = |0, 0.707, 0.707). The direction of the dashed lines in Fig.Q close to the tri- 
critical point corresponds to |Ai). The two dashed lines start at A* = 1, J|, one of them ends at the XY fixed point 
(A = 0, J 2 = 0) and the other ends in the Ising fixed point (A = 00, J 2 — 0) . Although these lines separate two 
different RG flows below and above the dash lines the flows represent a unique phase which is the spin-fluid for A < 1 
and the Neel phase for A > 1. These lines are supposed to be the Gaussian fixed points which have been reported in 
Ref.jjJ. In the spin-fluid and Neel phases, large distance physics is described by the Gaussian model where the system 
possesses the conformal symmetry |3j| . 

We have also calculated some critical exponents at the tri-critical point. In this respect, we have obtained the 
dynamical exponent, the exponent of order parameter and the diverging exponent of the correlation length. This 
corresponds to reaching the tri-critical point from the Neel phase by approaching A — > 1. The scaling of gap exponent 
is z ~ 0.83, the staggered magnetization goes to zero like Sm ~ |A — 1\@ with j3 ~ 0.82 and the correlation length 
diverges at A = 1 with exponent v ~ 2.27. The detail of a similar calculation can be found in Ref.[j|. 
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IV. SUMMARY AND DISCUSSIONS 

We have presented the second order RG approximation to obtain the phase diagram, staggered magnetisation 
and ground state energy of the XXZ chain with next-nearest neighbour interaction. We have observed that the 
second order RG equations must be considered to get a self similar Hamiltonian upon renormalisation processes. The 
ground state energy density of the nearest neighbour XXZ ( J2 — 0) model in the limit N — > 00 is then given by, 
e^ G = J2m=o 3^+1 ea(J^ m \ A^ m ^) where initially = J and A^ ) — A. We have considered the nearest neighbour 
case to be able to compare with the exact results. We have shown this quantity in Table. (1) where the results of 
the first and second order RG presented for < A < 0.9. The relative error has been derived by comparison with 
the exact result The second order correction improves the relative error at least to 65 percent. However the 

improvement of the 2nd order approximation is not monotonic for the whole range of A. 

An interesting result of the QRG scheme is the case of isotropic Heisenberg model (J2 = 0, A = 1) which runs to 
the frustrated next-nearest neighbour fixed point under the second order RG equations. This can be viewed as an ex- 
planation for the absence of long range order in the nearest neighbour antiferromagnetic Heisenberg model. Moreover, 
our primary calculation to take into account the higher order RG equations shows that the effective (renormalised) 
Hamiltonian will consist of the long range interactions. Taking into account the third order correction will result 
in an infinite range interaction of XXZ type. We would also like to mention a comment on the implementation of 
the 2nd order RG approach. The calculation is straightforward (not necessarily simple) whenever the zeroth order 
projection operator (Po) is composed of degenerate states. This makes possible to obtain the effective Hamiltonian 
in the analytic form. We have tried to get the 2nd order corrections for the S = \ XYZ model in the transverse field 
where the block ground-state is not degenerate. We were not able to obtain the analytic 2nd order corrections. This 
is important because the advantage of the QRG is the analytic form of results which helps to get the phase diagram 
even though it may not be quantitatively very accurate. If we are forced to do the RG approach in a numerical 
way then it is appropriate to implement the Density Matrix Renormalisation Group (DMRG) method. In that case, 
we loose the renormalisation group features of the running of the coupling constants but get the accurate numerical 
values. 
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